Towards jet tomography: 7-hadron correlations 
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Hard pQCD processes taking place in ultrarelativistic heavy-ion collisions are a well-calibrated 
probe. It is believed that the interaction with the surrounding medium of outgoing partons from 
a hard vertex is capable of revealing details of the medium. We demonstrate that correlation 
measurements of hard photon-hadron back-to-back coincidences are a tool suitable to extract such 
tomographic information. Introducing the concept of averaged evergy loss probability distributions, 
we first argue that almost no information about details of the medium evolution is reflected in the 
nuclear suppression factor Raa- Thus, a wide variety of jet quenching scenarios and geometries are 
compatible with the measured data. This problem can be overcome by a 7-hadron correlation mea- 
surement. We show that averaged probability distributions for quarks are accessible experimentally 
and sketch an analysis procedure capable of distinguishing different energy loss scenarios leading to 
the same nuclear suppression factor. 

PACS numbers: 25.75.-q,25.75.Gz 
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The expression 'jet tomography' often used to describe 
the analysis of hard pQCD processes taking place inside 
the soft matter created in an ultrarelativistic heavy-ion 
collision. In particular the focus is on the nuclear sup- 
pression of hard hadrons in A-A collisions compared with 
the scaled expectation from p-p collisions due to loss of 
energy from the hardparton by interactions with the soft 
medium (see e.g. [1 II II 0, If ) . 

The idea of jet tomography suggests that it is possible 
to make a 'cut' (greek to^lo^) through the outer layers of 
the medium and study directly the hot and dense core of 
the early evolution of the fireball created in the collision 
process. 

However, the nuclear suppression factor 
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is a rather integral quantity, arising in model calculations 
from the schematical convolution of the hard pQCD vac- 
uum cross section da AA c ~*f +x for the production of a 
parton /, the energy loss probability Pf(AE) given the 
vertex position and path through the medium and the 
vacuum fragmentation function D v f °^ h (z, fip) as 
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Here, fi/A(x, Q 2 ) denotes the distribution of parton i in- 
side the nucleus as a function of the parton momentum 
fraction x and the hard scale Q 2 of the scattering and 
&ij^f + k is the the partonic pQCD cross section 
As this expression has to be averaged over all possible 
vertex positions and paths for the emerging parton /, the 
question to what degree information about the structure 
of the medium is preserved and if such a measurement 
can be called tomography seems reasonable. 
In this paper, the question is investigated as follows: 
First, the concept of the averaged energy loss probability 
distribution (P(AE, E)) is introduced. We argue that 
this is a quantity which can be calculated in a wide vari- 
ety of models (and is thus suitable for model comparison) 
and is (for quarks) also experimentally accessible. We 
demonstrate that a wide range of different (P(AE, E)), 
corresponding to widely different physics mechanisms for 
jet quenching, yields a nuclear suppression factor com- 
patible with the measured data. Thus, it can be argued 
that there is little tomogrphic information in Raa- We 
show that the same variety of scenarios can easily be dis- 
tinguished by measuring the spectrum of hard leading 
away side hadrons correlated with a high px photon trig- 
ger and show to what degree tomographic information 
about the medium can be recovered from such a mea- 
surement. 



II. AVERAGED ENERGY LOSS PROBABILITY 
DISTRIBUTIONS 
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Considering the structure of Eq. one observes that all 
information about the medium is contained in the energy 
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loss probability Pf(AE) for a given parton / (while it is 
also possible to represent the energy loss of the leading 
hadron in the form of a medium-modified fragmentation 
function DJff h (z , fi 2 F ) , cf. e.g. it is more convenient 
to retain the representation for the present discus- 
sion). In a less schematic representation, the evaluation 
of Eq. @ hence requires spatial integrals over the pos- 
sible vertex positions, angular integrals over the direc- 
tion of propagation through the medium and line inte- 
grals along the path of the outgoing parton / through 
the medium (see e.g. @> for details). However, these 
spatial integrals factorize with the momentum-space in- 
gredients da^~* f+x and D v f "^ h (z, fj, F ). Thus, given a 
model for the medium and the energy loss, the integrals 
over all unobserved spatial quantities can be carried out 
separately, defining an averaged energy loss probability 
distribution (P(AE, E)). This quantity can in principle 
retain dependence on the initial hard parton energy E 
and contains all information about the medium convo- 
luted with the energy loss model. 

Such a quantity is also a highly intuitive way to char- 
acterize the medium: Given all assumptions about the 
medium and the energy loss, (P(AE, E)) immediately 
reveals which fraction of partons is completely absorbed, 
which fraction does not undergo energy loss at all and 
which fraction is shifted by some AE in the spectrum be- 
fore fragmentation occurs. In general, (P(AE, E)) will 
be probabilistic — either because the energy loss for a 
given parton path is treated on a probabilistic basis 
or, if this is neglected, because the distribution of in- 
medium pathlcngth needs to be taken into account and 
the initial vertex and path of an individual parton is not 
observed [iHTof. 

On the other hand, the averaged probability distribution 
for quarks is in principle accessible experimentally using 
the qg — > qj process HHH. If a photon with a large 
transverse momentum p7p is used as trigger, the momen- 
tum spectrum of correlated high p\ hadrons at angle « ir 
can be measured. The momentum spectrum can be cal- 
culated as 
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where we define the fraction of parton energy left after 
interaction with the medium as / = E ~ F E and trans- 
form (P(AE,E)) accordingly. For known D vac (z,Q 2 ), 
some information about the energy loss distribution can 
be inferred, dependent on the available momentum range 
and statistics. We postpone a more detailed analysis of 
this question to section Ivl 



III. SENSITIVITY OF Raa TO (P(AE,E)) 

Since, after doing the spatial averaging with the medium, 
(P(AE, E)) is the only place in Eq. (J2J where infor- 
mation about the medium enters, from this particular 
point of view a measurement of Raa only provides con- 
straints for the averaged energy loss probability distribu- 
tion. The degree of tomographical information in Raa is 
thus correlated with the sensitivity of Raa to details of 
(P(AE,E)). 

There are two main mechanisms which can reduce Raa- 
Absorption of a given fraction of partons ('downward 
shift' of the spectrum) and energy loss (without absorp- 
tion) of every parton ('sideward shift' of the spectrum). 
The true scenario will possibly contain any combination 
of the two. 

To illustrate this point, let us study some simple toy mod- 
els: 

• Geometrical suppression: We assume that the 
medium consists of an opaque core and a trans- 
parent halo. All partons in the core are absorbed 
and all partons in the halo escape without energy 
loss. In this case, 

(P(AE, E)) = TS(0) + (1 - T)S(E) 

(where T is the fraction of transparent region). 

• Smoothed geometrical suppression: Instead of tak- 
ing a sharp boundary, we allow for some transition 
region, i.e. we replace the (^-functions by Gaussians 
with a width d and replace the infinite opaqueness 
by some typical energy loss AEq Then, 

(P(AE,E)) = Ty^d/2exp[-E 2 /d 2 } 

+ (1 - T)yfirdexp[-(E - AE ) 2 /d 2 } 

• Typical energy loss: We assume that the whole 
medium is characterized by a typical energy loss 
AEq. In this case, 

(P(AE,E))=6(AE ) 

• Typical fractional energy loss: Here we assume that 
a hard parton typically loses a constant fraction 
(1 — /) of its energy. Then 

(P(AE,E)) = 5((l-f)E) 

• Extended semi-opaque medium: If the medium is 
large enough so that all partons are initially inside 
the medium and transparent enough such that sup- 
pression only accumulates for those partons which 
have a long in-medium pathlength, the suppression 
probability can be assumed to be flat as a function 



3 



0.3 



typical energy loss 
smoothed geometrical suppression 
semi-opaque medium 
hydrodynamical evolution (quarks) 
hydrodynamical evolution (gluons) 




0.75 



J- 0.5 



0.25 



10 

AE [GeV] 



• PHENIXdata 

typical fractional energy loss 

— typical energy loss 

— smoothed geometrical suppression 

— semi-opaque medium 

— hydrodynamics 
geometrical suppression 




P T [GeV/c] 



FIG. 1: Left panel: Averaged energy loss probability distributions for several toy models (see text, energy-dependent distribu- 
tions are not shown) and a full hydrodynamical simulation. Right: Nuclear suppression factor Raa for different toy models 
and the hydrodynamical simulation compared with PHENIX data. 



of AE up to some maximum value AE max (de- 
termined by the maximum in-medium pathlength) . 
We model this as 



(P(AE,E)) =iV/(l + exp 



AE - AE, n 



with N the normalization factor. 
Clearly, the physics (and the medium) underlying these 
scenarios is very different and any tomographic method 
should be able to resolve the differences. In Fig. ^ we 
show the different model probability distributions (in- 
cluding one extracted from a hydrodynamic model [l3j 
using the formalism outlined in All probability dis- 
tributions are shown with the convention that a parton 
is considered to be thermalized and absorbed into the 
medium as soon as E — AE = 0.5 GeV is realized. This 
determines the fraction of absorbed partons as a function 
of parton energy from the distributions. We also show the 
resulting Raa after folding with the parton spectrum and 
the fragmentation function. 

As apparent from Fig. ^ the calculated Raa is for most 
models well in agreement with the data beyond pt = 5 
GeV. In particular, the flatness of Raa is described well 
in most scenarios. The exception is the constant frac- 
tional energy loss where the calculation yield an Raa 
which continues to drop as a function of px- Thus, we 
may draw as a first conclusion that the flatness of Raa 
constrains probability distribution to have no strong de- 
pendence on the initial parton energy. The actual shape 
of the distribution is not constrained. Indeed the subset 
(P(AE)) performs rather well. 



IV. THE MAGNITUDE OF THE QUENCHING 

Naturally, an arbitrary probability distribution does not 
describe Raa — in each of the models, one (or at most 
two) parameters need to be adjusted. For the geometri- 
cal suppression, T = 0.18 is used (note that the resulting 



curve is not flat as naively expected as it still reflects 
the difference between nuclear and proton parton distri- 
bution functions) , the smoothed geometrical suppression 
requires T = 0.3 and AEo — 17 GeV (with a weak depen- 
dence of the result to d) , the typical energy loss requires 
AE = 5 GeV, for the typical fractional energy loss we 
find / = 0.7 and the extended semi-opaque medium leads 
to AE max — 14 GeV (with a weak dependence on d). 
What is common to all these parameter choices is that 
they make quenching in some sense substantial, either by 
setting a large scale for the typical shift of the spectrum 
or by making the transmission fraction T small. Thus, 
while the flatness of Raa favours energy independence, 
its magnitude indicates substantial quenching. 
Much effort has been made to condense this qualitative 
statement into a number, i.e. an averaged transport coef- 
ficient q which characterizes the typical opaqueness of the 
early medium. Different groups have reported widely dif- 
ferent results, e.g. 5-10 GeV 2 /fm (Wiedemann/Salgado), 
0.35-0.85 GeV 2 /fm (GLV) , « 2 GeV 2 /fm (AMY) and 
3-4 GeV 2 /fm (Majumder) 0). Based on these re- 
sults, (dis)-agreement with pQCD expectations has been 
claimed. 

The main difference in the extraction of q in the frame- 
works cited above seems to be the question how to define 
the averaging procedure for a medium at which the trans- 
port coefficient is in principle different at each point in 
space and time and only a fraction of the partons may 
emerge from the medium while others are absorbed (and 
hence can't carry any information), or even to define a 
typical average pathlength of a parton in the medium. 
The same problem is apparent from a different angle in 
Fig. ^ — why should one expect to be able to extract a 
typical quenching power of the medium when widely dif- 
ferent energy loss probability distributions (correspond- 
ing to very different maximal and typical quenching and 
very different average pathlengths) result in the same 
Raa? 

For this reason, we have in [8j adopted a different proce- 
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dure: Taking_o = 2 • K ■ e 3 / 4 with K = 1 as the pQCD 
expectation we can for a given model of the space- 
time dynamics fit K to the measured Raa and investigate 
the disagreement with pQCD by doing the full spatial 
integration in Eq. J3J) instead of resorting to schematic 
averaging procedures. We found that, within a class of 
evolution models which describe bulk hadronic observ- 
ables such as WT-spectra for 7r, K, and p and dN c h/dr] 
(which are hence substantially more constrained than 
typical models used to calculate Raa) K can still vary 
from close to 1 to about 5, dependent on detailed assump- 
tions about the longitudinal and transverse flow profile 
and the magnitude of a s . 

Thus, the conclusion is that even the magnitude of Raa 
does not allow detailed conclusions about the early evo- 
lution dynamics. This, along with the independence of 
Raa to the detailed shape of the averaged energy loss 
probability, indicates that there is almost no tomograph- 
ical information to be gained from Raa alone and that 
more differential observables need to be studied in order 
to learn about the early evolution. 



V. PHOTON-HADRON CORRELATIONS 

The idea underlying the use of 7-hadron correlations to 
distinguish between different scenarios is very intuitive: 
If we had a monochromatic parton source and could mea- 
sure the spectrum of outgoing partons, the effects of ab- 
sorption of a fraction of partons and a sideward shift due 
to some energy loss would be completely different. In 
the case of absorption, the intensity of the source would 
be reduced but it would still be monochromatic, in the 
case of a shift the intensity would be unchanged but the 
spectrum would be altered and in general no longer be 
monochromatic. A photon trigger can provide in prin- 
ciple the necessary monochromatic parton source and 
a study of the momentum spectrum of hard away side 
hadrons carries the momentum information. Since the 
photon escapes from the fireball without further interac- 
tion, there is no geometrical trigger bias and the distri- 
bution of hard vertices for observed triggers is just given 
by the nuclear overlap as for Raa- 

The main practical problem is that the outgoing hadron 
fragments and we have to solve Eq. to gain access 
to the desired energy loss probability distribution. With 
the left hand side known at a finite number of points, 
with finite statistics and in a limited momentum range, 
this is not an easy task. 

We study the problem by embedding gq — > jq events at a 
given photon energy of 15 GeV into the Monte-Carlo sim- 
ulation designed to study dihadron correlations |l6j |. By 
replacing the dynamical energy loss calculation with the 
toy model probability distributions described in section 
IIIII (this is possible as the vertex distribution of events 
is the same as for Raa, see above) we can simulate the 
momentum spectra of hard hadrons back to back with 
the photon for each of the energy loss models. 



In order to extract the energy loss mechanism, we try to 
solve Eq. (@J with a trial ansatz. In the case of geomet- 
rical suppression, we can write 



(P((f,E 7 )) = T-8(l) + (l-T)-6(0) 



Inserting this into Eq. (@J yields 
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thus, the momentum spectrum divided by the fragmen- 
tation function is expected to be flat if geometrical sup- 
pression is realized. We show this representation of the 
simulated data for the different scenarios in Fig. [21 (we 
do not discuss constant fractional energy loss further, as 
this can be already identified on the level of Raa)- 
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FIG. 2: Momentum spectrum of hard hadrons correlated 
back-to-back with a photon trigger, normalized to the expec- 
tation of geometrical absorption (see text). 

The individual models can clearly be told apart. To 
go beyond that and recover some degree of information 
about the underlying probability distribution, we make a 
trial ansatz where we have in addition to the transmittion 
('T') term a typical energy loss ('L') 



(P((f,Ej) = T-6(l)+L-S((pl-AE )/pl) + (l-T-L)-6(0) 



^7 dN I r.vac/Ji /„7 ' 



Using the normalized quantity p^-^-/D vac (pl^/p^), our 
trial ansatz yields 



„7 dN 



D< 



D' 



m 



T + L- 



vac ( Pt \ 
\Pt J 



(8) 



Clearly, a simultaneous fit of T, AEq and L is not a 
promising strategy as the first two terms in Eq. (jTJ may 
become degenerate for AE — > 0. Thus, we chose to ad- 
just T to the lowest data point in each sample and fit L 
and A_Bo accordingly 
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The resulting parameters characterize the individual sce- 
narios more or less well. For the geometrical suppres- 
sion we find T = 0.19, I ~ and hence AEo undefined 
as expected. For the apparently rather similar hydro- 
dynamics calculation, we find T = 0.12, L = 0.07 and 
AEo = 4.38, i.e. in addition to ~ 12% partons without 
energy loss there are some ~ 7% partons undergoing on 
average about 4 GeV energy loss observed, the rest is 
absorbed. 

Let us next discuss the typical energy loss scenario. Here 
we find T — and L = 1.53, AE = 5.24. While the en- 
ergy loss value put into the simulation is approximately 
recovered, the result for L is a surprise as we would ex- 
pect T, L < 1 and the sum of transmission, loss and ab- 
sorption term to be unity. However, due to the finite 
momentum range considered here, we do not see the full 
distribution and contributions at low momentum can bal- 
ance the conservation of probability. 
For the smoothed geometrical suppression the fit yields 
T = 0.2, L = 0.25 and AE = 4.28. This is rather similar 
to the results of the semi-opaque medium where we find 
T = 0.16, L = 0.29 and AE Q = 3.6. 

Thus, while the ratio of L/T provides some measure to 
characterize the different scenarios as absorption or loss- 
dominated, it is instructive that the average value of en- 
ergy loss 'seen' in this study is always close to 4 GeV. 
This is presumably rather related to the assumed cuts 
- for a 15 GeV parton, significantly more energy loss 
does not result in much hadronic yield above 4 GeV, sig- 
nificantly less does not result in the correct amount of 
quenching unless there is an absorption term. Thus, by 
varying the spread between trigger and associate hadron 
momentum one opens the possibility to probe more and 
more of the energy loss probability distribution. 

VI. SUMMARY 

We have argued that averaged energy loss probability dis- 
tributions are a useful tool to compare individual model 



calculations, as they characterize an energy loss scenario 
according to the degree of transmission, absorption and 
energy loss of partons. Making several toy models for dif- 
ferent pictures of energy loss, we have demonstrated that 
Raa is incapable of constraining details of the energy 
loss probability distribution; a more differential observ- 
able is needed instead. Even the average magnitude of 
the quenching observed is not a useful means to charac- 
terize the actual conditions in the medium as it is un- 
clear how the averaging procedure should be defined in 
a model-independent way. Unfortunately, Raa is, due 
to its insensitivity to details, unsuitable to provide guid- 
ance. 

As a possible means to overcome this problem and to 
access the averaged energy loss probability distributions 
we have suggested 7-hadron correlations. We have car- 
ried out a proof of concept study and argue that such a 
measurement is able to make a distinction between the 
different toy models studied here and also to character- 
ize them to some degree. The same information is in 
principle accessible via dihadron correlations, cf. [lfij . 
However, in this case the analysis is much complicated 
by the fact that the initial energy of the partons, given a 
high pt trigger hadron, is only known on a probabilistic 
basis and by the additional problem that energy loss on 
the trigger leads to a systematic distortion of the ver- 
tices leading to a triggered event. While the latter is 
desireable from the point of view of studying pathlength 
dependence, it is clearly not advantageous for studying 
the momentum distribution of the associated hadron in 
a clean way |l7| . 
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